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Effects of strongly selective additives on volume phase transition in gels 
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We investigate volume phase transition in gels immersed in mixture solvents, on the basis of a three-component 
Flory-Rehner theory. When the selectivity of the minority solvent component to the polymer network is 
strong, the gel tends to shrink with an increasing concentration of the additive, regardless of whether it is 
good or poor. This behavior originates from the difference of the additive concentration between inside and 
outside the gel. We also found the gap of the gel volume at the transition point can be controlled by adding 
the strongly selective solutes. By dissolving a strongly poor additive, for instance, the discontinuous volume 
phase transition can be extinguished. Furthermore, we observed that another volume phase trasition occurs 
far from the original transition point. These behaviors can be well explained by a simplified theory neglecting 
the non-linearity of the additive concentration. 



I. INTRODUCTION 

Swelling behavior of a polymer network was first in- 
vestigated by Flory and Refiner— . On the basis of this 
work, discontinuous volume phase transition (VPT) of 
gels was predicted theoretically 2 . Thence, the VPT has 
received much attention from both scientific and indus- 
trial viewpoints 3 — . In various types of gels, addition of 
solutes into the solvent often plays very important roles 
in their volume change 3 - - — . For example, the VPT was 
first realized in a polyacrylamide (PAA) gel immersed in 
a mixture of water and acetone 3 -. In this experiment, the 
solvent quality, or x parameter, is controlled by changing 
the volume fraction of the acetone. It is also well known 
that addition of salts dramatically affects the behaviors 
of ionic gels 7 - - — . Furthermore, a nonioninc hydrated gel 
such as poly(A^-isopropylacrylamide) gel (NIPA) changes 
its volume discontinuously. In the hydrated gels, the 
VPT is attributed to temperature-dependence of hydro- 
gen bondings of the polymer network^S. It has been 
reported that its VPT is affected by adding salts^i - — 
and other additives^ - —. 

It is often assumed that the composition of a mixture 
solvent is constant in and out of a gel. This assump- 
tion is referred to as single liquid approximation (SLA). 
The SLA is experimentally confirmed in a NIPA gel in a 
dimethyl sulfoxide-water mixture, for instance^ 2 .. How- 
ever, the compositions of mixture solvents are not nec- 
essarily homogeneous 2 ^ - —. There, the difference of the 
composition between the interior and exterior of the gel 
would lead to a dramatic effect on its volume change. 
Some extended Flory-Huggins models dealing with a 
polymer network and two spieces of solvent molecules 
were developed to investigate the volume change of a 
polymer network in a mixture solvent 2 ^ - —. In these the- 
ories, the volume change is characterized by the compo- 
sitions of the three components and three x parameters. 
They were studied systematically for fixed sets of the x 
parameters. Usually, the x parameters in polymer sys- 
tems depend on environment parameters such as tem- 
perature. Since it is quite complicated to incorporate 
the temperature dependences of the x parameters, the 



swelling behavior in mixture solvents has not been fully 
understood. 

Recently, the effects of solutes with strong selective 
solvation on phase behaviors of a water-oil mixture was 
investigated 2 ^ - —. By dissolving a very small amount of 
strongly hydrophobic solute into the water-oil mixture, 
the oil-rich phase containing the solute is precipitated 
even in the one-phase region. This finding implies that 
such solutes would affect the VPT in gels. 

The aim of this study is to clarify the effect of the ad- 
ditives on the VPT in gels, based on a simple theoretical 
argument. In particular, we focus on the cases, in which 
one of the solvent components has strong selectivity to 
the polymer network. 



II. MODEL FREE ENERGY 

Swelling behaviors of gels are well described by the 
Flory-Rehner theoryi. The free energy F consists of a 
mixing part F m [ x and an elastic part F c \ as 



F„ 



(1) 



For a gel immersed in a mixture solvent, the mixing free 
energy F mix is given by 2 ^ -25 ' 27 ' 29 



-Fmix = k B Tv Q 1 [V g f(<f>l g , 02g, 03 ) + V s f((j> ls , 2s , 0)] , 

(2) 

where fee is the Boltzmann constant and T is the tem- 
perature. V g and V s are the volumes of the interior and 
exterior of the gel, respectively. In this model, we as- 
sume that the volumes of each monomer of the polymer 
network and the solvent molecules are equal to a char- 
acteristic volume vq for simplicity. /(0i,02,03) is the 
Flory-Huggins type mixing free energy of the ternary sys- 
tem as 3 -! 
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Here, 0i and 02 are the volume fractions of the first and 
second components of the mixture solvent, respectively. 
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In Eq. ([2]), 0^ g and 0i S stand for the volume fraction of 
the i-th component (i — 1,2) in and out of the gel. 03 
is the volume fraction of the polymer network. Since 
all the polymer chains are chemically connected form- 
ing a single network, the polymers are not dissolved out- 
side the gel, that is, 03 S = 0. In Eq. ([3]), the first two 
terms stem from the translational entropy of the solvent 
molecules. Here, the translational entropy of the net- 
work is neglected. The last term is the interaction en- 
ergy among the solvent molecules and the monomers of 
the network. Xij is the interaction parameter between 
the i-th and j-th components. 

We employ Flory's rubber elasticity as22, 
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where V^o is the volume of the gel in a reference state, 
v is the density of cross linking points, and 7^ is the 
elongation ratio in the A:-axis (k = x,y, z). B is a nonlin- 
ear elastic coefficient, which is often assumed B = 1/2. 
In polyelectrolyte gels, the translational entropy of the 
counter-ions is renormalized into B as B + b, where b 
is the number of the dissociable monomers per effective 
chain. In this study, however, we do not consider the elec- 
tric charges explicitly. In the case of isotropic swelling, 
the elongation ratio is coupled with the volume fraction 
of the polymer network as 
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where 030 is the volume fraction of the network in the 
reference state. Therefore, the elastic energy Eq. (jH) is 
rewritten by a function of 03 as, 
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We impose the incompressible conditions for both in and 
out of the gel, 



!»lg + 02g + 03 = 1, 



(7) 
(8) 



01s + 02s = 1- 

Here, we define a grand potential as 

= F miK + F cl - Li2{V g (/>2 g + V B <t> 2s ) + k(V s + V B ), (9) 

where /i2 and k are Lagrange multipliers to conserve the 
amount of the second component and the total volume. 
The equilibrium state is characterized by minimizing i7, 
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Then, we obtained the equilibrium conditions: 



At(</>2g,</>3) = A(02s,O) 
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where jl is the reduced chemical potential given by 
A(02>3) = In-; -7— — — 

1-02-03 
+ Xl2(l - 202 - 03) + (X23 - X3l)03- (14) 

Eq. (fT2^) and ([13"]) represent the balances of the chemi- 
cal potential for the second component and the osmotic 
pressure, respectively. 

Ho is a part of the osmotic pressure for the gel without 
the second solvent component. It stems from the mixing 
free energy and given by 

n o = -03-X3i03-m(l-03)- (15) 
fl e i is the contribution of the elasticity as 




If the solvent does not contain the additive 02 = 0, the 
balance of the osmotic pressure is expressed by n -l-II el = 
0. _ 

n a( j is the contribution of the second component and 
is given by 
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where G = X23 — X31 — X12 is a parameter describing the 
affinity of the second component to the polymer network. 
n and n a( jd come from the mixing free energey F m i x . If 
G < and \G\ ^> 1, the additive tends to be adsorbed 
selectively to the polymer network. If G > 1, on the 
other hand, the additive would be expelled from the gel. 

In this study, we assume that the binary solvent is 
completely mixed outside the gel. We note that the vol- 
ume of the equilibirated gel does not change if we add 
more mixture solvent whose composition is equal to the 
equilibrated outer solvent. Hence, we take a limit of 
V s /V g — > 00 with fixing 02 S , so that 02 S is uniquely deter- 
mined for a fixed /i2- Hereafter, we use 02 S as a control 
parameter. 



III. RESULTS AND DISCUSSIONS 
A. The first volume phase transition 

We numerically study the effects of additives on vol- 
ume changes of gels. First, we focus on the cases, in 
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FIG. 1. The swelling curves of a gel immersed in a mixture 
solvent. We set 3O = 5.0 x 10~ 2 , vv = 1.0 x 10~ 2 , B = 
0.75, and X12 — 0.0. In the absence of additives, the gel 
undergoes volume phase transition with increasing X3i- We 
dissolve solutes of X23 = —1.0 (pro-gel) in (a) and X23 = 2.0 
(anti-gel) in (b). The black solid lines are calculated from 
Eq. 42T). 



which the gel can undergo VPT without addtives. We 
set 3O = 5.0 x 10~ 2 , vvv = 1.0 x 10~ 2 and B = 0.752. 

We vary X31 continuously to induce the volume change 
with fixing the other parameters. We obtain the swellng 
curves of the gel by solving Eqs. (TT21 and (TIB")) numeri- 
cally. Figure 1 shows the swelling curves for several val- 
ues of 02s- In the absence of additives, the VPT occurs 
at X31 = X3it(~ 0.93). In Fig. 1 (a), we dissolve the ad- 
ditive of X12 = 0.0 and X23 = —1-0 , which has affinity 
to the polymer network (pro-gel). As its concentration 
02s is increased, the transition point shifts to higher X3i- 
In Fig. 1 (b), the additive of X12 = 0.0 and X23 = 2.0 
is dissolved into the solvent. This additive dislikes the 
polymer network (anti-gel). It is shown that the transi- 
tion point shifts to lower X31 with an increasing 2s - For 
the cases as in Fig. 1, the SLA can explain the effects 
of the solute on the volume changes qualitatively. In the 
SLA, the volume fractions of the solvent components in- 
side the gel are assumed to be 0i g = (1 — 03) (1 — 02s) 
and 02 g = (1 — 03) 02s- Then, the interaction parameter 
between the polymer network and the mixture solvent is 
approximated as 



X31 



X3l(l - 02s) + X2302s - Xl202s(l ~ 02s) 

G0 2s + O(0L). (18) 
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This means that the interaction parameter changes effec- 
tively with 02s depending on G. For the plotted range 
in Fig. 1 (a) (x3i S [0.92,0.95]), G remains negative, 
so that the solvent becomes more good, swelling the gel 
with an increasing 02 S - As shown in Fig. 1 (b), on the 
other hand, the solvent changes to more poor and the gel 
shrinks with 02 S for the solute of positive G. 

In Fig. 1, the absolute value of the resulting G is rather 
small (|G| < 2). Next, we study the effects of solutes of 
strong selectivity \G\ ^> 1. Figures 2 (a) and (b) show the 
swelling curves in the mixture solvents of X23 = —9.0 and 
X23 = 11-0 , respectively. The other parameters are the 
same as those in Fig. 1. In the both cases, the transition 
points are shifted to lower X31, regardless of whether the 



FIG. 2. The swelling curves of a gel immersed in a mixture 
solvent. We set 3O = 5.0 x 10~ 2 , vv = 1.0 x 10~ 2 , B = 
0.75, and X12 = 0.0. We dissolve solutes of X23 = —9.0 (pro- 
gel) in (a) and X23 = H-0 (anti-gel) in (b). The black solid 
lines are calculated from Eq. (|21[) . The additives have strong 
selectivities to the polymer network. 
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FIG. 3. (a) The shifts of the trantision point from xiit 
are shown for several values of X23- The shift is numerically 
obtained with the same parameters as in Figs. 1 and 2. For 
X23 = 9.0 and 11.0, the transition points are terminated since 
the discontinous volume changes disappear, (b) The depen- 
dence of Ax3it on X23 is plotted for 02s = 0.01. The broken 
curve represents g(cf>3 C ,G) [Eq. ((26)) ] with 03 C = 0.239 and 

X31I = 0.93. 



solute is good or poor. This behavior is in contrast to the 
volume changes in Fig. 1 and indicates the SLA does not 
work well when the additives have the strong selectivity 
to the polymer network. 

Figure 3 (a) shows the dependences of the transition 
shift Ax3it on the additive concentration 02 S , where 
AX3U = XSH - and *£> is the transition point 
without additives. An increase of Ax3it is observed for 
X23 ~ —1-0, whereas Ax3it is lowered with an increas- 
ing 02s for the other additives. In the plotted range of 
02s, Ax3it has linear dependence on 02 S . The transition 
shifts for 02s = 0.01 are plotted with respect to X23 in 
Fig. 3 (b). It indicates a non-monotonic behavior of the 
transition point. 

The difference of the additive effects between X23 = 
— 1.0 and X23 = —9.0 is notable, since the both addivies 
have tendencies to be adsorbed onto the polymer net- 
work. When |x23 1 is small, the difference between 02 g 
and 02s is rather small. As explained by Eq. (TT51) . the 
addition of the solute with negative X23 makes the mix- 
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FIG. 4. The gap of the volume fraction 03 at the transition 
is plotted with 2s . We set 3 o = 5.0 x 10~ 2 , uvo = 1.0 x 10~ 2 , 
B = 0.75, and X12 = 0.0 as in Figs. 1 and 2. 

ture solvent more good to the polymer. Thus, the gel 
is simply swollen with an increasing 2s - If |%23 1 is large 
enough, on the other hand, the mixture solvent inside the 
gel becomes much denser than that outside the gel. By 
shrinking its volume, the polymer network tends to in- 
crease the contact points to the additive molecules. This 
non-linear effect would give rise to the difference between 
X23 = -1-0 and X23 = -9-0. 

Figure 2 also shows that the gap of the volume tran- 
sition depends on the additive concentration when the 
additive has strong selectivity. As 02 S is increased, the 
volume gap is increased for X23 = —9.0 [see Fig. 2 (a)], 
while it is decreased for \23 — 11-0 [Fig. 2 (b)]. For 
X23 = 11.0, in particular, the gap disappears eventually 
at a certain concentration 02s = 02st , above which the gel 
does not undergo the VPT [see below] . Figure 4 plots the 
gap of the volume fraction A03t(= 03t+ — fat-) at the 
transition point, where 03t+ and 03t- are the volume 
fractions just below and above the transition point. For 
a large positive value of X23 , Afa t decreases toward zero 
with an increasing 02 S - 

It is known that the elasticity parameter B influences 
the nature of the VPTl. Fi gure 5 shows the swelling 
curves of a gel of B = 0.72 with varying X3i- In the ab- 
sence of additive, the volume of this gel changes contin- 
uously. By dissolving the pro-gel additive of X23 = —9.0, 
the swelling curve becomes non-monotonic with respect 
to X31, so that the gel undergoes the VPT. Figures 2 (b) 
and 5 indicate that we can arbitrarily induce or erase the 
VPT by adding solutes with the strong selectivities. 

We assume that only a subtle amount of the additive 
is dissolved in the outer solvent, i.e., 02s <C 1. From 
Eq. (fl"2|) . the volume fraction of the additive inside the 
gel is expressed by 

fa g = (l-fa)e- G ^fa s + 0(4> 2 2s ). (19) 
Substituting Eq. (JT5J) into Eq. (fT7|). we obtain 

n ad = ln(l - 2s ) - ln(l - fa s e~ G ^) 

+Gfa(l-fa)e~ G ^fa s + 0(tf 2g ) 

= [e- G ^{\ + Gfa(l - 3 )} - 1] 02s + 0(4>l g )- 

(20) 
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FIG. 5. The swelling curves of a gel immersed in a mix- 
ture solvent. We set fao = 5.0 x 10~ 2 , uvo = 1.0 x 10~ 2 , 
B = 0.72, and X12 = 0.0. We dissolve the pro-gel solute of 
X23 = —9.0. The black solid lines are calculated from Eq. (|2ip 
. When 02s = 0, the gel changes its volume continuously with- 
out VPT. 



Then, the swelling curve is approximatedly given from 
Eqs. (H5j), (HSJ and flU by 

X3i = -0 3 ~ 2 ln(l - 3 ) - afa 5/3 + {fi- 1)0 3 - 1 

+ [e- G ^{l + Gfa(l - 3 )} - l]03~ 2 ^2s, (21) 

1 /3 

where a — vvo/<j> 30 and (3 — vv^BjfaQ. Since the right 
hand side of Eq. (f2l"j) includes X31 via G, we cannot obtain 
an analytical solution of X3i- Then, we replace G to 
G = X23 — xiit — X21 and regard G as a fixed parameter. 
The approximated curves of Eq. (|2 1[) are drawn in Figs. 1, 
2 and 5. They are well in agreement with the numerical 
solutions. 

Regarding /? as a variable, a tri-critical point of the 
VPT is given by 

5X31 ^ Q d 2 X3i =Q (22) 

^ 03c, /3c 03c, /3c 

Here, a is treated to be fixed. By solving Eq. (|22|). we 
obtain 03 C and /3 C , which give the swelling curve passing 
through the tri-critical point. In the vicinity of the tri- 
critical point, the swelling curve is approximated by 

X31 ~ X31t - —72— ~ (fa - 03c) + 7^(03 - 03c) 3 , (23) 
03c 6 

where u is a positive constant. If j3 < f3 c , the curve 
°f X3i changes monotonically with 03, such that the gel 
volume changes continuously with changing X3i- On the 
other hand, inflection points appears when /3 > j3 c . If 
so, the gel becomes mechanically unstable, so that it 
exhibits the discontinuous volume changes as A0 3t = 
2yJ(P - f3 c )/u/fa c at xsi = XsiA 

The part of the osmotic pressure ll a( j shifts the tri- 
critical point in Eq. (|23l as 

X3it=X^t+5(03c,G)02s, (24) 
P c = pW+h(fa c ,G)<f>l c fa s . (25) 
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FIG. 6. The contours of the functions g [Eq. ([26j] an d h [Eq. 
(|27|l ] are drawn in G-03 planes (a) and (b), respectively. The 
transition point is increased linearly with 02s in the region of 
for positive g. In the region of negative h, the volume gap of 
the transition point is increased. 
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FIG. 7. The swelling curves of a gel immersed in a mixture 
solvent in a wider range of X3i- We set 030 = 5.0 x 10 -2 , 
vv = 1.0 x 10~ 2 , B = 0.75, and xia = 0.0. We dissolve the 
strongly good solute of X23 = —16.0 and 02s = 2.0 x 10~ 7 . 
The black solid lines are calculated by Eq. (|21|) • The second 
volume phase transition is induced. 



Here, fi c gives the tri-critical point in the case of 02s = 
0. From Eq. (f^Tj) . the prefactors g and h are given by 

5 (03, G) = 03 V G03 U + G0 3 (1 ~ 0s)} - 1], (26) 
G) = 3 " 3 [2 - e-^iG^Kl - 3 ) + 2G0 3 + 2}] . 

(27) 

g and h change their signs depending on 03 and G. Their 
dependences are shown in Fig. 6. If g((f>3,G) is positive, 
the transition point shifts to higher xai with an increas- 
ing 02s and vice versa. Eq. (1261) for 02s = 0.01 is drawn in 
the broken curve in Fig. 3 (b). Here, we set 03 C = 0.2392 
and X3it = 0.93. This approximated curve is in ag- 
greement with the numerical solution. When 03 C <C 1, 
g m — G(G + 2)/2 and it changes its sign at G = and 
G = — 2. A large negative h enhances the discontinuity of 
the volume change as observed in Fig. 2 (b). If 03 C <C 1, 
h is expressed as h « G 2 (G + 3)/3 and it changes its 
sign at G = —3. Thus, we can possibly induce the VPT 
even in a gel, which originaly shows a continuous vol- 
ume change, by adding solute of — G ^> 3. The gap 
of the VPT disappears as A03 t oc V02st — 02s, where 
02st = {P — /^c }/i _1 03~ c 2 . This disappearing behavior is 
observed in Fig. 4. 

In this paper, only the numerical solutions for X12 — 
0.0 are presented. However, we confirmed that the essen- 
tially same features are observed for any set of Xl2 and 
Xai if the resultant G is the same. 

B. The second volume phase transition 

In Fig. 7, we plot the swelling curves in a wider range of 
X3i- We set 030 = 5.0 xl0~ 2 , uv = 1.0 xl0~ 2 , B = 0.75, 
X12 = 0.0 and X23 = —16.0. It is indicated that another 
discontinuous volume phase transition can occur at X31 
larger than %3it corresponding to the first transition. We 
confirmed that this second phase transition is observed 
even if v — (or, a = ft = 0), while the first one disap- 
pears. This fact indicates that this second instability is 



independent of the network elasticity and has a physical 
mechanism different from those for well-studied volume 
phase transitions 3 . 

Scott reported that various types of phase diagrams are 
realized for ternary mixtures (polymer solutions in binary 
mixtures)^. They can have several critical points, below 
which three phases coexist. Analogous to the ternary 
mixtures, we consider that the second volume transition 
observed in Fig. 7 is attributed to the bulk instability of 
the mixing free energy. Since there are some differences 
between ternary mixtures and gels in mixture solvents, 
the Scott's argument cannot be simply applied to the 
gel systems. The most important differences are that 
a gel has the elasticity and never exhibits a one-phase 
homogeneous state. 

We have not obtained a simple explanation on the sec- 
ond transition, since non-linearities, which are hard to 
treat analytically, would play imporant roles in the bulk 
instability. In Fig. 7, we also draw the swelling curve 
obtained by Eq. (f2"Tj) . which neglects the non-linearity of 
02s- The approximated curve also exhibits the second 
transition, although it does not coincide quantitatively 
with the numerical solutions for large 03. However, this 
fact implies that the non-linearity of 02 has a minor con- 
tribution to the essential mechanism of the second tran- 
sition. 

As discussed above, the gel becomes mechanically un- 
stable when 9x3i /903 < 0. In Fig. 8, we plot dx^i /dfa 

and — /i(03, G)02s, where X3?(^3) is the swelling curve of 
the gel in the solvent of the first component. In the range 

satisfying 9xii /903 < — /i(03, G)02 S , the gel would ex- 
hibit a discontinuous volume change. It is indicated that 
strong non-linearity of h((f>a,G) with respect to ^3 is a 
possible origin of the second transition. 
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FIG. 8. The theoretical curves of dxsi /d(j>3 and 

—h(4>3,G)(f>2 B are shown. When dxfi/dfo < —h((j> 3 ,G)(j)2 B , 
the gel becomes mechanically unstable, showing the second 
volume phase transition. 



IV. SUMMARY AND REMARKS 

We studied the volume phase transition of gels 
immersed in binary mixtures, based on the three- 
component Flory-Rehner model. Assuming that the vol- 
ume fraction of the second solvent component is small, 
we reformulated the Flory-Rehner model into a simple 
model with a new parameter G{= x?3~ X3i~Xi2)- From 
numerical solutions and the simplified theory, we found 
the following behaviors of the volume phase trasition of 
the gel. 

(i) When the selectivity of the second component to the 
polymer network is small, the composition of the mixture 
solvent in the gel is close to that out of the gel. Here, 
the single liquid approximation works well. The renor- 
malized interaction parameter between the solvent and 
the network depends linearly on the composition of the 
binary mixture. Thus, we can swell the gel by dissolv- 
ing the additive good to the network, whereas, the gel is 
shrunken when the poor additive is dissolved. 

On the other hand, the difference of the composition 
becomes large when the additive selectivity is strong. 
Owing to the non-lineariry of the composition difference, 
the gel tends to be shrunken as the additive concentra- 
tion is increased, regardless of whether it is good or poor. 
From the simplified theory, the dependence of the tran- 
sition point on the addtive concentration is found to be 
proportinal to c*X3it/<902s ~ 9{4>3c, G), which is approxi- 
mated as g « —G{G + 2)/2 for a hyper swollen gel. 

(ii) By dissolving the strongly poor additive, we can ex- 
tinguish the discontinous volume phase transition. Fur- 
thermore, we can induce the volume phase transition in 
the gel, which does not exhbit the discontinous gap in the 
absence of the additive, by dissolving the strongly good 
additive. In the vicinity of the (first) tri-critical point, 
this behavior is well described by the simplified theory 
as shown in Fig. 6 (b). 

(iii) Far from the transition point, we found that an- 
other volume phase transition can occur. This second 
instability is observed even if the elasticity is negligible. 



The transition is caused by the mixing instability. Our 
theory indicates that the non-linearity of <fe , which is well 
described by Eq. (|2T|) . does have a major contribution to 
the transition. 

Here, we make some remarks about our results. 

1) If only van der Waals interaction is taken into ac- 
count, the interaction parameters should be positive, i.e., 
Xij > 0. However, a variety of molecular forces, such as 
ion-dipole interaction, hydrogen bonding and hydropho- 
bic interaction, would also influence the volume phase 
transition of the gel. When studying their contributions 
explicitly, we have to consider the microscopic degrees 
of the freedom in a more specific manne r 21 ' 33 i 34 . In this 
work, we assume that the macroscpic \ parameter can 
have a large negative value, by renormalizing these mi- 
croscopic degrees into them. In a future work, we would 
like to investigate the connections of such microscopic 
interactions and our phenomenological theory. Effects of 
the selective salts 2 ^ on polyelectrolyte gels are also inter- 
esting. 

2) There are a lot of experiments dealing with NIPA 
gels and PAA gels to study the effects of the additive on 
the volume phase transition. Most of additives shrink 
NIPA gels in agreement with our results. Since only a 
few experiments have reported the difference of the ad- 
ditive concentrations between in and out of the gel s 22 ' 35 , 
however, we cannot verify the quantitative validity of our 
model yet. Further experimental studies are highly de- 
sired. 
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